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SONIC RADIATION FROM A CIRCULAR PISTON 
WITH IMPULSE EXCITATION 


By Leon Bess 

Electronics Research Center 


SUMMARY 

Closed analytic expressions have been derived for the sonic 
pressure as a function of space and time when a circular piston 
is excited so that its velocity varies with time as a Dirac delta 
function. These expressions are valid not only for the far field 
but also for the near field. Moreover a number of typical plots 
of the sonic pressure versus time for various special points have 
been made using the results of the calculations. The basis of 
the method of calculation is the approximation of the circular 
piston by an array of very small spheres uniformly distributed 
over the area of the piston. Each of these spheres are assumed 
to have their radius increase as a step function in time at the 
instant of excitation. The sonic field for such an excited 
sphere is well known and is relatively simple in form. 

******* 

The sonic field for a circular disk with impulse excitation 
has been solved for distances large compared to the disk radius 
(ref. 1). For small distances and distances comparable to the 
radius, the evaluation becomes too complex when conventional 
methods are employed. 

In the work to be reported here, a different approach to the 
solution has been taken which allows its evaluation for both the 
near and far field without too much difficulty. The first step 
in attaining the solution is to illustrate the configuration of 
the geometric system and this is done in Figure 1. As can be 
seen, the flat excitor disk lies in the XY plane with its center 
at the origin. The field measurement point, M, lies in the YZ 
plane, and has the coordinates (y,z). It is assumed that all 
points in the front face of the disk move in the + z direction 
during excitation, and that these points move with a velocity, u, 
given by: 


u = A6 ( t) (la) 

As can be seen from studying Eq. (la) , the whole disk face 
suddenly moves a distance. A, in the z direction at t=0, and then 
remains stationary thereafter. (6(t), incidentally, is the Dirac 
delta function) . 






For the method adopted here, the starting point of the whole 
calculation is the determination of the radiation field from a 
very small sphere. The radius of the sphere is to suddenly in- 
crease by an amount, Ag. The assumption is now made that a very 
large number of these small spheres with their centers distri- 
buted uniformly over the area of the disk shown in Figure 1 would 
approximate the action of a disk whose thickness suddenly in- 
creases by the amount, 2A, at t=0. One disk face moves in the 
+z direction by an amount. A, and the other in the -z direction 
by the same amount. The justification of this assumption forms 
the basis for the validity of the whole treatment and will now 
be gone into in detail. 

The desired expression for the sonic pressure, P 3 , for the 
present disk problem must satisfy only two conditions. The first 
is that it be a solution of the differential equation usually 
referred to as the wave equation (i.e., - ( 1 /c^) 

P^ = 0). The second condition is that P^ satisfy the following 
relations in the time interval, 0<t<e(2/C). 

P^ = C 6 (t ± z/c) for y <_ a (lb) 

P^ = 0 for y > a (lc) 

C is a constant, a is the disk radius, C is the velocity of sound, 
and e is to be chosen so that it is a very small number (typically 
having a value around 1 / 100 ). 

The radiation field from a small sphere has been treated in 
several places (ref. 2). (This result is only valid if the radius, 
r Q , of the small sphere is such that A<g<_r 0 <<cTp where Tp is the 
duration of the excitation pulse.) This condition can be satis- 
fied in most practical cases. The result for the sonic pressure, 
Pg, as a function of space and time is: 

P S ( ^' t} = dt [ 6(t “ r/G) J (ld) 

where y is the density of the fluid, and S is the surface area of 
the sphere. 

Since the sonic pressure, Pg, for a single small sphere is a 
solution to the wave equation everywhere except on the disk, the 
total sonic pressure from a large array of these spheres will be 
such a solution also, because it is merely the linear super-posi- 
tion of a large number of the functions, Pg, (and the wave equa- 
tion is a linear differential equation). Thus, the total array 
pressure satisfies the first condition mentioned above. It 
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satisfies the second condition also, but this cannot become 
apparent until the total array pressure is calculated. The cal- 
culation has actually been performed in the treatment that follows. 
By temporarily skipping ahead, it can be verified that the total 
array pressure (to be labeled *D> satisfied the boundary condi- 
tions at times near t=0 . Thus, it can be seen that the P D , as 
determined from Eqs . (9b) and (19) show that it agrees with Eq. 

(lc) . Thus, since P D satisfies the two conditions of acceptance 
given above, it follows that P D = P^, the desired solution. 

From the description of the method of calculation given 
above, it would appear that the sonic pressure, P^, generated by 
the action of the disk can be given by the following expression: 


2tt 


p D (r,t) 


= (q/S) 


If 


p o(r,t ; p,cj>) pd<t>dp 


( 2 ) 


2 2 1/2 

where p = [x + y ] ; tan 4> = (x/y) ; and q is a source strength 

constant to be determined. 


In evaluating the integral of Eq. (2) , the (p integration is 
to be performed first. The result of the cp integration is dp D 
which is the pressure wave that would result if only a ring of 
radius, p, and width, dp, were excited. Using Eqs. (2) and 
(Id), the explicit form of dp D is: 




dt 


2v 


I 


(dcp/r) 6 (t - r/c) 


where r = 


+ y2 + p2 _ 2yp cos (f) 


1 1/2 


(3) 


The evaluation of Eq. (3) involves using the properties of 
a Dirac delta function. It can thus be seen that dp D /0 only for 
those values of t where t = r/c. Moreover, since only r appears 
explicitly in the argument of the delta function and the variable 
of integration is <f>, a factor of c | dr/d<j> | ~1 must appear. The 
result of the cp integration of Eq. (3) is then: 


dP = 0 for t < x and t>T 


dP T 


ry A s p dp\ d 


4tt 


dt l r | dr/dcf> 




(4a) 

(4b) 
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A .m 2 2 7 ' T = , r~2 2 2 7 

where x =Vz + y + p + 2yp v Vz + y + p - 2yp 


V 2 2 2 

z + y + p - 2yp cos cf> t 

Note that there are two values of <j) that can satisfy the (j). 
equation. From Eqs. (3) and (4b) it follows that: 

|r~| = py | sin <J> t | = j V4p 2 y 2 - (D 2 - p 2 ) 2 (5) 

d> = <(. t 

, 2 _ 2,2 2 2 

where D =ct -z - y 

Let r) be defined so that: 

2 2 2 2 2 „ , - , 
c t = z + y + p - 2pyn (6) 

It then follows from Eqs. (4b) and (5) that: 



It is to be noted that there is an extra factor of 2 in 
Eq. (7) to account for the fact that there are two values of <j> 
satisfying the <f> t equation. n is a function of t only with 
ri = -1 when t = x; rt = 1 when t - x and - l<ri<l when x<t<x. 
Equation (7), therefore, provides a more suggestive form for dP D . 
Finally, it is of interest to note that: 



With the use of Eqs. (7) and (8) it is possible to see that 
for large distances from the ring (i.e., where z 2 + y 2 >> p 2 ) the 
time and space dependence of dP^ begins to resemble that of Eq. 

(Id), the pressure wave from a small sphere. This is what is to 
be expected and therefore provides an accuracy check on the 
calculations . 
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The radiation pressure, P D , for the total disk can now be 
obtained by integrating the result for the ring (i.e. dP D given 
by Eqs. (4a), (4b) and (5)) with respect to the variable, p. When 

this is done, the result can be presented in the form: 

P = 0 when t>t and t<t (9a) 

1 J V* 



for t < t < t (9b) 

The limits of integration, p and p will in general depend on t as 
well as y. The exact determination of these limits is rather 
complicated and, therefore, must be deferred until later. The 
times, t and t, are also to be determined. 

It is possible to perform the p integration in Eq. (9b) by 
changing variables from p to a so that a = - p2. if the 

quantity in the square bracket in Eq. (9b) is designated as I, 
then it can be shown that after substituting o for p, I has the 
form: 


I = 




do 

2 2 
4y a - a 


( 10 ) 


The integral of Eq. (10) can be evaluated exactly and the result 
is : 


I 


arc sin 

2 2 2 
D + 2y - p 

•»/ 2 2 

( 

2y Vd + y 


1 

2 


arc sm 


D 2 + 2Y 2 - p 2 


2y Vd 2 + y 2 


(ID 


There now remains the task of evaluating the limits, p and p. 
It can be seen that in doing this, two distinct cases naturally 
arise. Case I is where y lies in the range, 0<y<a, and Case II 
is where y>a. Case I is to be considered first and this can best 
be done by referring to Figure 2. This is seen as a representa- 
tion on the XY plane of the measurement point, M, and the disk 
trace (shown by the cross-hatched bar) . There is also a repre- 
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sentation of the disk on the YZ plane. Upon studying this figure, 
it is apparent that the point on the disk closest to the point, M, 
is the point, L, which is the intersection of the line, LM, with 
the surface of the disk. The line, LM, lies in the XY plane and 
is parallel to the x axis. It is obvious that the minimum time, 
t, in the pulse is determined by the product of (1/c) and the line 
length, LM. Similarly, it can be seen that the point on the disk 
having the greatest distance to M is N, the point at the bottom 
of the disk. Thus, the maximum pulse time, t, is (1/c) times the 
length, MN. From the above considerations it follows that: 

t = z/c (12a) 

t = (1/c) Vz 2 + (y + a) 2 (12b) 


It can also be seen from a study of Figure 2 that for any 
given t (in the range t to t) , the locus of all points on the 
surface of the disk having a distance, ct, to the point, M, is a 
circle whose center is at the point, (0,y,0) and which has a 
radius E, given by 



z 


2 * 


It follows from this that the minimum radius, p, of any ring con- 
tributing to the integral of Eq. (9b) is determined by the short- 
est distance from the origin to any point on the locus circle. 
Similarly, the maximum radius, p, is determined by the greatest 
distance of any point on the locus circle to the origin. 


From these considerations 

and a study of 

Figure 2 , 

it can be 

shown that p and p are 

V 

P = 

given by 

y + E 

for E < a - y 


(13a) 

A 

P = 

a 

for E > a - y 


(13b) 

P = 

y - E 

for E < y 


(14a) 

P = 

? - y 

for E > y 


(14b) 

2 

P = 

2 2 
D Z + 2y 

- 2 y Vd 2 + y 2 

for all E 

(14c) 


where E = Vc 2 t 2 - z 2 = Vd 2 + y 2 


Using the results of Eq. sets (13) and (14) in Eq. (11) , the 
final result for I when y<a is: 
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I = ir/2 


I = 


for t < t < t 
2 


IT 1 

j - 2 arc sm 


„ 2 2 
+ 2y - a 


2y 


Vd 2 


+ y 


(15a) 

(15b) 


where t = (1/c) Vz 2 + (a - y) 2 for t < t < t 

For Case II where y>a, it is possible to determine p and j3 
in the same manner that they were determined in Case I. v The main 
difference is that instead of using Figure 2 to determine the 
various geometric relations. Figure 3 must now be used. Just as 
in Case I, the length, MN, represents the longest distance between 
any disk point and the point, M. However, for Case II, QN repre- 
sents the shortest distance instead of LM. From an analysis of 
Case II similar to the one already performed for Case I, it is 
possible to show that the following relations are valid: 


t = (1/c) Vz 2 + (y - a) 2 

s/ 


(16a) 

t = (1/c) Vz 2 + (y + a) 2 


(16b) 

P = y - 5 for E, < y 

V 


(17a) 

P = £ - y for ? > y 


(17b) 

p 2 = D 2 + 2y 2 - 2y Vy 2 + D 2 

for all 5 

(17c) 


where 


€ = Vd 2 + y' 


p = a 


for t < t < t 


(18) 


Using the results of Eq. sets (17) and (18) to obtain p and 
the value of the quantity, I, of Eq. (11) can be seen to be (for 
y>a) : 


I = ( j - j arc sin 


2 2 2 
D + 2y - a 


2y Vd 2 


4- y 


(19) 


for t < t < t 


With the derivation of Eq. sets (15) and (19), the quantity, 
I, has been completely determined and, therefore, with the use of 
Eq. set (9), the form of the pressure, P D , can be derived. To 
completely determine P , however, it is necessary to find the 
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Figure 3.- Case 


y > a 


CASE I ; Y < a 



Figure 4.- Plot of I vs t for Cases I and II 





value of the source strength constant, q. To do this and also to 
provide a check for some of the calculations presented above, a 
simple special case of the sonic radiation field is now to be 
treated. This case is where y = 0 (i.e. M lies on the X axis) . 
Here it is obvious that 


t = z/c and t 
Moreover it follows that: 


1 

c 



+ 




( 20 ) 


If the variable of integration is changed from p to r, and 
the relation pdp = rdr is used, Eg. (20) becomes transformed to: 


P 


D 


2tt 



for t < t < t 

V 


( 21 ) 


The integration indicated in Eq. (21) can be readily per- 
formed and it produces a pair of step funtions . These become 
delta functions again when differentiated with respect to t so 
that the final result for P D is: 



r - 


( / n o \ 

P D = 

(1/2) qcyA g 


6 (t - z/c) — 6 ft — (1/c) \z z + a J 

/\ _ 


for t < t < t 


The result of Eq. (22) is now to be subjected to two quali- 
tative checks. It is to be expected that when a-*- 00 , the form of P 
approximates that where the piston is an infinite plane with all 
the surface points moving with the velocity defined by Eq. (la) . 
The solution for this can readily be obtained and is found to have 
the form of the delta function, 6 (t - z/c). It can be seen that 
as a-*- 00 , the form of Pq does indeed approach the required form of 
6 (t - z/c), since the second term in Eq. (22) appears only after 
a very long time. Thus, Eq. (22) is able to pass the first check. 


As a second check it is to be expected that when z-»-°°, P D 
approaches the form obtained from a pulsating spherical monopole 
as given by Eq. (Id). It is to be noted, first, that 
[ 6 (u+A ) -6 (u) =Ad/du [ S (u) ] ] when A<<u. Thus it follows from 
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Eq. (22) when z>>a: 


P 


D 


§^ A s 


2 

a 

2 z 


d 

dt 


6 t 


z_ 

c 


(23) 


Since, in this case z=r, it can be seen that Eq. (23) agrees in 
form with Eq. (Id) and, thus Eq. (22) is able to pass the second 
check also. 


In order to determine q, the source strength constant, the 
same boundry conditions used for the small sphere (ref. 2) can 
be invoked here. These are that on the surface of the piston the 
relation, y du z /dt = -dP^/dz, is valid. Using Eq. (la) to give 
u„ and Eq. (22) to give P D , it follows that the boundry condi- 
tions require that the various parameters be related as follows: 



(24) 


Since the general expression for Pq given by Eq. (9b) involves 
the entire quantity, (qA s ). the source strength constant, q, is, 
in effect, determined by Eq. (24) . 

It is of interest to note that the form of Pp given by 
Eq. (22) agrees with the Pq given in the main text (from Eqs. 

(9b) , (15a) , and (15b) ) when y = 0. It can be shown that for 

these forms to agree, I in Eq A set (15) must-be a function such 
that 1 = 0 for t < t and t > t; I = tt/ 2 for t > t > t. This con- 
dition is satisfied by Eq. set (15) since it_readily ^follows from 
Eq. (12b) and Eq. set (15) that when y = 0, t = t. Thus, Eq. (22) 
is able to pass the final check. 

The final task of the treatment is to present the quantity, 
dl/dt, in explicit form. Doing this requires a straight-forward 
but rather tedious differentiation of the quantity, I. Without 
presenting any of the calculational details, the final result is 
given by the following equations: 

= 0 for t > t and t < t (25a) 

dt v 


dl _ 1 
dt 2 


for t < t < t 


ft «=- 


J) 


Vl^ 


(G - J) 


(25b) 
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where : 


G = 


2,2 2 , 2 

_ c t - z + y 


2y^ 22 2 


fc*t‘ - z‘ 


Case I; y < a. 


dG _ c 
dt 


» I 2.2 2 2 ] 

tLct - z -yj 

r 2~2 2 \ 3/2 

c t - z J 


2y (c 


J = 


„ x f2~2 2 

2y vc t - z 


for t < t < t 


(26a) 


dJ 

dt 


2 2 . 
-act 


_ , 2,2 2 3/2 

2y \c t - z ' 


for t < t < t 


(26b) 


Case II; y > a. 


J 


2y Vc 2 t 2 



dj. = 
dt 


2 2 

-act 


2y 


2,2 

c t 


- z 2 W 2 


(27a) 


(27b) 


With the derivation of the quantity, dl/dt, the sonic radia- 
tion field has now been completely determined. The character of 
the field can probably be best obtained by substituting in the 
pertinent equations above to get a plot of I versus t for differ- 
ent values of x and z. This operation has been carried out, and 
the results are given below. However, it is felt to be more 
instructive to first obtain a rough qualitative description of 
the sonic field from the study of the pertinent equations. 


Before proceeding with the description, it is necessary that, 
as an aid to the discussion, a rough general plot of I versus t 
for both Cases I and II be made. This is done in Figure 4. It 
can be seen from this figure that it ijs not necessary to calculate 
dl/dt for the time interval, t < t < t, in Case I, since it always 
has the same form in this interval. 


From Eq. (22) it can be seen that when y = 0, Pq is in the 
form of two delta functions desplaced in time with respect to one 
another. This is to say that when the piston is excited by a 
very sharp time pulse, the pressure waves it generates on the x 
axis will also be sharp pulses. Off the x axis but for y < a. 
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the P D versus t plot remains two pulses displaced in time and the 
leading pulse stays sharp (i.e. it is a delta function). The 
lagging pulse, however, gets spread out (to a pulse duration of 
around a/c) . This result is perhaps the most important one in 
the whole treatment since it accounts for the fact that sharp 
pulse excitations can produce sharp sonic pulses over extended 
regions . 

When the measuring point, M, is far off the x axis (where 
y > a) it is found that both the leading and lagging pulse is 
spread out with pulse durations of the order of ac. Moreover, 
as y becomes larger, the amplitude of the pulse (as well as that 
of the quantity, I) decreases in value. 

As mentioned above, the derivation obtained here has been 
used to prepare plots of I versus t and (dl/dt) versus t for 
various spatial points. Figures 5, 6, and 7 give I versus t and 
Figures 8, 9, 10, 11, and 12 give (dl/dt) versus t. These curves 
are self explanitory. Figure 12 essentially give? a picture of 
the "wave front" at the instant of time, t = (4.0) (a/c). The 
vertical arrows in the above figures indicate Dirac delta 
functions . 



Figure 5.- I vs t 
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(c/a) t = 4.0 



Figure 12.- dl/dt vs t (t = (4.0) (a/c) 


17 



REFERENCES 


1. Morse, P. M. , and Ingard, K. U. : Theoretical Acoustics. 
McCraw-Hill Book Co., New York, 1968, pp. 337-338. 

2. Ibid: pp. 309-310. 


18 


NASA-Langley, 1971 23 C“140 



National Aeronautics and Space Administration 
Washington, D. C. 20546 




OFFICIAL business 

PENALTY FOR PRIVATE USE $300 


FIRST CLASS MAIL 



POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AN 
i ADMINISTRATION 


06U 001 48 
AIR FORCE WEA 
KIRTLAND AFBf 


ci 3 ns 71088 00903 

PONS laboratory /wlol/ 

NEW MEXICO 87117 


ATT E. 


LOU BOWMAN, 


CHIEF, TECH. LIBRARY 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



